Phase ordering with a global conservation law: Ostwald ripening and coalescence 
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Globally conserved phase ordering dynamics is investigated in systems with short range corre- 
ct . lations at t = 0. A Ginzburg-Landau equation with a global conservation law is employed as the 
phase field model. The conditions are found under which the sharp-interface limit of this equation 
is reducible to the area-preserving motion by curvature. Numerical simulations show that, for both 
critical and off-critical quench, the equal time pair correlation function exhibits dynamic scaling, 
and the characteristic coarsening length obeys l(t) ~ t 1 ^ 2 . For the critical quench, our results are 
in excellent agreement with earlier results. For off-critical quench (Ostwald ripening) we investigate 
the dynamics of the size distribution function of the minority phase domains. The simulations show 
that, at large times, this distribution function has a self-similar form with growth exponent 1/2. The 
O ' scaled distribution, however, strongly differs from the classical Wagner distribution. We attribute 
this difference to coalescence of domains. A new theory of Ostwald ripening is developed that takes 
into account binary coalescence events. The theoretical scaled distribution function agrees well with 
that obtained in the simulations. 
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I. INTRODUCTION 



Phase ordering is emergence of order from disorder through domain growth and coarsening. The standard setting 
when phase ordering occurs is a temperature quench from a high-temperature disordered phase into a two-phase or 
a multiphase region. Phase ordering has been the subject of extensive research during the last two decades fl. An 
important simplifying assumption in phase ordering theory is dynamic scale invariance. According to this assumption, 
the coarsening system possesses, at late times, a single relevant dynamic length scale lit) (the characteristic domain 
size) which grows with time as lit) ~ t a |Q. It is by now well established that in systems with short range correlations 
a = 1/2 for non-conserved (model A) dynamics, while a = 1/3 for locally conserved (model B) dynamics . 

There is, however, an important additional coarsening mechanism: globally conserved phase ordering [^|, |3|, ^, [|, [Fj. 
\ Globally conserved dynamics can be thought of as model A dynamics constrained by global conservation of the 
jy-j ■ order parameter: for example, Ising model with fixed magnetization. This global conservation law is maintained 
y—{ by an external field (for example, a magnetic field) which depends on time but is uniform in space. The globally- 
■ conserved phase ordering is accessible in experiment. Consider the sublimation/deposition dynamics of a solid and 
its vapor in a small closed vessel kept at a constant temperature below the melting point. As the acoustic time 
scale in the gas phase is short compared to the coarsening time, the gas pressure (and, consequently, density) remain 
uniform in space, changing only in time. This character of mass transport in the vapor phase makes the coarsening 
dynamics conserved globally rather than locally. An important characteristics of globally-conserved dynamics is 
interface-controlled kinetics, in contrast to the bulk-diffusion-controllcd kinetics typical for locally-conserved systems. 
Interface-controlled kinetics was investigated in the context of growth of small platinum particles supported on alumina 
substrates in an oxidizing environment There are additional examples of cluster growth on surfaces |j, where 
it was found possible to single out the interface-controlled kinetics [|l0|. There is also strong evidence in favor of 
O ■ globally conserved interface-controlled transport during the coarsening of clusters in granular powders driven by a 
^ ■ low- frequency electric field jyj |i~2| . 

, Part of the theoretical importance of globally conserved phase ordering lies in the fact that it enables an access to 
off-critical quenches in the simpler model A dynamics (with global conservation). Thus, it allows one to determine 
which characteristics of the system depend on the volume (or area) fraction e and which do not. 

Dynamic renormalization group arguments show that global conservation should not change the growth law [ jl3| . 
This early result was confirmed by particle simulations with short range correlations in the initial conditions: for 
critical (e = 1/2) p[ fT4[ and for off-critical (e < 1/2) Q quench. Recent phase field simulations of systems with 
long-range (power-lawjcorrelations in the initial conditions have also shown dynamic scale invariance with the same 
growth exponent a = 1/2 ||. Therefore, a = 1/2 independently of e. On the other hand, the autocorrelation function 
P, |l5| and persistence exponent jL5| were found to be e-dependent. 

Globally conserved dynamics are related to a wide range of multiphase coarsening systems. Sire and Majumdar 0] 
showed that in the large-g limit the dynamics of the g-state Potts model are equivalent to the dynamics of the globally 
conserved model with an area fraction e — 1/q. The large-q limit of the Potts model is of practical importance as it 
describes correctly some of the dynamic characteristics of dry soap froths |16| and of the coarsening of polycrystalline 
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materials |L7fl . 

In the limit of a vanishing volume fraction of the minority phase the late stage of coarsening is describable by the 
mean field theories of Ostwald ripening. Lifshitz and Slyozov |1| developed such a theory for the bulk-diffusion- 
controlled (or locally conserved) dynamics. They showed that the size distribution function of the minority domains 
approaches, at large times, a self-similar form. Correspondingly, the average size of the minority phase domains grows 
with time like i 1 / 3 . Following the seminal work by Lifshitz and Slyozov [|l8[, Wagner developed a similar mean- field 
theory for the interface-controlled (globally conserved) Ostwald ripening |19|. The Wagner's theory yields a growth 
law t 1 / 2 for the average domain size, and a different (broader) shape of the scaled distribution function. We shall refer 
to this scaled distribution function as the Wagner distribution. More recently, it was shown that interface-controlled 
Ostwald ripening appears in the sharp-interface limit of scalar Ginzburg-Landau equations (and its modifications) 
with a global conservation law (2| || [|, || p0| . 

Although the simple theories of Lifshitz-Slyozov and Wagner were developed more than 40 years ago, they are 
still very useful in phase ordering theory. For example, Sire and Majumdar [Q employed the Wagner distribution to 
calculate the equal-time pair correlation function. Lee and Rutenberg |ll| used the two theories of Ostwald ripening 
for calculating the autocorrelation exponent and persistence exponent for the locally- and globally-conserved systems. 

Many works were devoted to extensions of the Lifshitz-Slyozov theory to finite volume fractions. Already Lifshitz 
and Slyozov [[l8) made an attempt to go beyond their simple model and account for coalescence. Later it became 
clear that, in the locally-conserved systems, the dominant effect unaccounted for by the simple theory is inter-domain 
correlations, rather than coalescence. At small area fractions, the relative role of correlations is of order e 1 / 2 
while the relative role of coalescence is of order e. Therefore, an account of coalescence without a proper account of 
correlations is an excess of accuracy. 

The situation is quite different in globally-conserved systems, and this fact has not been recognized until now. 
Correlations between neighboring domains are exponentially small in this case ||. Therefore, coalescence is expected 
to give the dominant correction to the theory of Wagner |T^]. We shall report numerical simulations that show a 
strong effect of coalescence at moderate e. Specifically, we find that, at large times, the size distribution function of the 
minority domains has a self-similar form with the "normal" growth exponent 1/2. The scaled distribution function, 
however, strongly differs from the Wagner distribution. We attribute this difference to coalescence and develop a new 
theory of Ostwald ripening that takes coalescence into account. 

The outline of the rest of the paper is the following. In Section II we briefly review the phase field model of globally 
conserved phase ordering: a scalar Ginzburg-Landau equation with a global conservation law. The sharp-interface 
asymptotic limit of this equation is introduced and reduced, in 2D, to a simpler model of area-preserving motion 
by curvature. The criteria for the validity of this reduction are obtained and presented in Appendix. The model of 
area-preserving motion by curvature is used to obtain dynamic scaling laws for the characteristic coarsening length 
and for the effective magnetic field. The results of the phase field simulations for the critical quench (e = 0.50) and 
off-critical quench (e = 0.25) are presented and analyzed in Section III. In Section [V a new theory of Ostwald ripening 
is developed. The theory leads to a non-linear integro-differential equation for the scaled distribution function. The 
solution of this equation is in good agreement with the scaled size distribution function found in the phase field 
simulations. In Section ^ we summarize our results. 



II. PHASE FIELD MODEL AND SHARP-INTERFACE LIMIT 



Globally conserved phase ordering dynamics are describable by a simple phase field model 
the free energy functional has the Ginzburg-Landau form: 



F[u] 



2 fi 



(Vm) 2 + V(u) + Hu 



d d r. 



In this model 



(1) 



and the dynamics follow a simple gradient descent: 
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d t u = —\i —— = 8 V 2 u + a (u - 
ou 



H(t)) . 



(2) 



Either no-flux or periodic boundary conditions can be used. In Eqs. (|l]) and (||) u(r, t) is the coarse-grained order 
parameter field, V(u) — (1/4) (1 — u 2 ) 2 is a symmetric double- well potential, 8 is the diffusion coefficient, ^ is the 
characteristic rate of relaxation of the field u to its stable equilibrium values, and d is the dimension of space. The 
effective uniform magnetic field H (t) changes in time so as to impose the global conservation law: 



(u{r, t)) = L~ d J u(r, t) d d r = const , 



(3) 
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where L is the system size and the integration is carried out over the entire system. Integrating both sides of Eq. 
over the entire system and using Eq. (J3|) and the boundary conditions we obtain: 



Therefore, Eq. (0) takes the form: 



dtu — 5 V 2 u + /i (u — u 3 ) — pb{u — u 3 ) 



(4) 



(5) 



a globally-constrained Ginzburg-Landau equation (GLE). From now on we shall concentrate on the 2D case. 

At late stages of the coarsening process, the system consists of domains of "phase 1" (where u is close to —1) and 
"phase 2" (where u is close to 1) separated by domain walls. The domain walls can be treated as sharp interfaces 
H H|, as their characteristic width A = (5/n) 1 / 2 is much smaller than the characteristic domain size l(t) which grows 
with time. At this stage H(t) is already small, H(t) <C 1, and slowly varying in time. The phase field in the phases 
1 and 2 is almost uniform and rapidly adjusting to the value of H(t), so u ~ — 1 — H(t)/2 and u ~ 1 — H(t)/2 in 
the phases 1 and 2, respectively. Under these conditions, the so called "sharp- interface theory" holds. The normal 
velocity of the interface v n is given by [|| : 



v n (s,t) = S K (s,t) + (6fi) 1 / 2 gH(t), 



(G) 



where s is a coordinate along the interface, n(s,t) is the local curvature, and g = — 3/V2- A positive v n corresponds 
to the interface moving towards phase 1, while a positive k corresponds to an interface which is convex towards phase 
2. The dynamics of H(t) are described by ||: 



H(t) 



4A(t) 
L 2 



Sn(s,t) + (5n) 1/2 gH{t) 



Here n(s,t) is the interface curvature averaged over the whole interface: 

1 



k(s, t) 



A(t) 



k(s, t) ds , 



(7) 



(8) 



and A(t) = <f> ds is the total perimeter of the interface. Equations 
formulation for the GLE with a global conservation law. 
Let us denote by A(t) the total area of phase 2: 



and (|7|) provide a general sharp-interface 



A(t) 



d 2 r. 



tt(r,t)>0 



Equations (|J) and (]?]) can be used to calculate the rate of change of A(t): 

A{t) = I v n (s, t)ds = A(f) \Sn(s,t) + (<5/i) 1/2 gH(t) 



Using Eqs. ([?]) and ^ we obtain: H = AA/L 2 which yields the global conservation law: 



A(t) 



L 2 H(t) 



— const . 



(9) 



(10) 



The second term in Eq. (|10|) corresponds to the bulk order parameter u being biased by H(t). One can use Eq. ( |10| ) 
instead of Eq. (0) in the general sharp-interface formulation of the problem. 

In some important cases this formulation can be simplified further || ^, |2^] . When the two terms on the right hand 
side of Eq. (0) approximately balance each other, 



H(t) 



Ifi 

9 \M 



1/2 



K,(s,t) , 



(11) 



the area of each of the two phases remains constant. In this case Eq. ([jj) takes the form: 



V n ( s ,t) = S K(s,t) — K,(s,t) 



(12) 
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Dynamics ( ]l2| ) are known as area-preserving motion by curvature in 2D, and as volume-preserving motion by mean 
curvature in 3D §, H |||. Due to the presence of the non-local term k this model is different from the Allen-Cahn 
equation 25 v n = 6 k, which represents the sharp-interface limit for non-conserved (model A) dynamics 

A simple example where the area-preserving dynamics cannot be used is the dynamics of a single circular domain 
of the minority phase in a "sea" of the majority phase || |2(J. Another example is the dynamics of a "donut": a 
single domain of the minority phase with an inclusion of a majority phase domain [ ^6| . Therefore, the first question 
we need to address concerns the general conditions under which the area- preserving dynamics, Eq. (|l2]), represent an 
accurate approximation to the more general sharp-interface theory, Eqs. (|6|) and (|7]). These conditions are derived in 
Appendix. 

Now we employ the area-preserving model and do simple dynamic scaling analysis. (In the rest of the paper we 
are using dimensionless variables and put 5 = fJ, = 1.) For critical quench we have H(t) = and k(s,t) = because 
of symmetry between the two phases (we neglect finite-size effects). Therefore, the globally conserved dynamics for 
critical quench are identical to the non-conserved (model A) dynamics. Using the Allen-Cahn equation v n {s,t) = 
n(s,t), one arrives at the well-known scaling law l(t) ~ t 1 / 2 Q. 

Turning to the off-critical quench, we notice that, under the scaling assumption, the interface velocity can be 
estimated as v n ~ all /alt. Each of the two terms on the right hand side of Eq. (12|) is of order l/l(t). Equating and 
integrating, we again obtain l(t) ~ t 1 / 2 . Therefore, global conservation does not change the dynamic scaling for any 
area fraction. This result was previously obtained by dynamic renormalization group arguments applied to Eq. (|5|) 
(with a Gaussian white noise term) 113], and by particle simulations B. For the off-critical dynamics of H{t) we have: 
\H(t)\ = | < K(s,t) > /g\ ~ l/l(t) ~ t- 1 ' 2 . 

Though the dynamic exponent is independent of the area fraction, other characteristics can depend on it. In the 
following Section we report numerical simulations that address area-fraction-dependent quantities. 



III. NUMERICAL SIMULATIONS 

We performed extensive simulations by directly solving Eq. (|J) with initial conditions in the form of "white noise". 
The simulations were done for two different values of the area fraction of the minority phase: e = 0.50, and e = 0.25, 
corresponding to a critical and off-critical quench, respectively. In both cases the results were averaged over 10 different 
samples. Eq. was discretized and solved on a 1024 x 1024 domain, with mesh size Ax = Ay — 1 and periodic 
boundary conditions. The coarsening process was followed up to a time t = 3000. An explicit Euler integration 
scheme was used to advance the solution in time, and the Laplace operator was discretized by second order central 
differences. A time step At = 0.1 was required for numerical stability. The accuracy of the numerical scheme was 
monitored by checking the (approximate) conservation law ([l0]) of the general sharp-interface theory. It was found 
that this conservation law is obeyed with an accuracy better than 0.02% for t > 4 and better than 0.008% for t > 30 in 
the critical quench case. In the off-critical quench the approximate conservation law ( |l0| ) is obeyed with an accuracy 
better than 1% for t > 30. To avoid any misunderstanding here and in the following we notice that, in all cases, the 
integrated order parameter [see Eq. (^)] is conserved exactly by the numerical scheme. 

It is convenient to introduce an auxiliary density field p(r, t) = (1/2) [u(r, t) + 1]. The minority phase is identified as 
the locus where p(r, t) > 1/2. Typical snapshots of the coarsening process are shown in Fig. |l| for the critical quench, 
and in Fig. |^ for the off-critical quench. For the critical quench the system consists of interpenetrating domains of the 
two phases. For the off-critical quench the morphology is that of Ostwald ripening ^?J: larger domains of the minority 
phase grow at the expense of smaller ones. As the minority phase area fraction is not very small, binary (and even 
triple) coalescence events are clearly seen in Fig. ||. Overall, the coarsening morphologies resemble those observed for 
locally-conserved system: in numerical solutions of the Cahn-Hilliard equation ps| and in particle simulations p9| . 
An important difference is an apparent absence of correlations between neighboring domains in Fig. ^[ 

To analyze the coarsening dynamics, the following quantities were sampled and averaged over the 10 initial condi- 
tions: 

1. The area of phase 2. 

2. The circularly averaged equal-time pair correlation function: 

< P (r>,t) P (ri + r,t)>-< P (ri,t)>> 
<p 2 (r',i) > -<p(r',t) > 2 

3. The characteristic coarsening length scale l(t), determined from the condition C(l,t) = 1/2. 

4. The effective magnetic field H(t) computed from Eq. (H). 
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5. The size distribution function of the minority phase domains (for the off-critical quench). 

For critical quench we found that the area of the minority phase is constant with an accuracy better than 0.03% 
at all times. The situation is quite different for the off-critical quench. Here there is a systematic trend in the area 
fraction of the minority phase. Still, with time this quantity approaches a constant value. Deviations from this 
constant value become less than 3% for t > 100. This approximate area conservation plays a crucial role in the theory 
of Ostwald ripening (see Sec. ^v|) . 

Figure ^ shows, on a single graph, the scaling forms of the correlation function C(x), where x = r/l(t), for the 
critical and off-critical quench. The /(t)-dependence is presented in Fig. |J. A comparison of the scaling forms C(x) 
with those obtained in particle simulations of globally conserved |6j and non-conserved pofl dynamics for critical 
quench is also shown. The three curves for the critical quench almost coincide. For off-critical quench, the C{x) 
curve is slightly different from the curves for the critical quench. A similar weak dependence of the scaled correlation 
function on the area fraction of the minority phase was observed in locally conserved systems [^8[ |3lJ . 

Figure [| shows corrected power-law fits l[t) = Iq + bt a which yield a — 0.50, Iq — 0.5, and b = 1.2 for the critical 
quench, and a — 0.51, lo — 1.3, and b = 0.9 for the off-critical quench. A pure t x l 2 power-law line serves as a 
reference for the expected late-time dynamic behavior. Therefore, l(t) obeys the expected t 1 ^ 2 dynamical scaling 
law, in agreement with the predictions of the dynamic renormalization group analysis p3[ and area-preserving sharp- 
interface theory. The difference in the values of the amplitudes b again indicates a dependence on the minority phase 
area fraction. 

The time history of l/\H(t)\ for the off-critical quench is presented in Fig. |5|. The data is fitted by a corrected 
power-law: l/\H(t)\ = a + ct a with a = 7.4, c = 2.3, and a = 0.51. Also shown is a 2.3 i 1 / 2 power law, serving as a 
reference to the expected late-time dynamics. We conclude that |i?(t)| ~ t" 1 / 2 , as predicted by the sharp-interface 
theory. The significance of the value of the amplitude c will be discussed in Section [iy| in the context of our theory of 
Ostwald ripening with coalescence. One can distinguish in Fig. || small "fluctuations" of l/\H(t)\ around a smooth 
trend. (This is in contrast to the /(t)-dependence where no fluctuations are observed.) To interpret these fluctuations 



we use Eqs. ( |l l|) and ( J50[ ) to obtain: l/\H(t)\ ~ A(f)/iVa(t). A(t) is a continuous function of t, whereas A^i) behaves 
discontinuously at the time moments when domains disappear due to shrinking and merging events. Thus H(t) serves 
as a "domain counter" . 

For critical quench, H (i) exhibits very small irregular fluctuations around zero. The typical values of H (t) in this 
case are of the order of 10 -5 , and we interpret these fluctuations as finite-size effects. 

As we have shown, the scaled correlation function only weakly depends on the area fraction. A much more sensitive 
diagnostics of the off-critical quench dynamics is provided by the size distribution function of the minority phase 
domains. We found that, at late times, this function exhibits dynamic scaling. Figure ^| shows the scaled form $ num 
of the distribution function obtained in the simulations with the GLE. The scaled variable on the horizontal axis of 
Fig. ^ is £ = R/t 1 / 2 , where the effective radius of each domain is defined as R — {Ad/'n) 1 ^ 2 and Ad is the domain 
area. Function $ lmm was obtained, at each moment of time, by multiplying the values of the distribution function, 
found in the simulations, by i 3 / 2 . The dynamic exponents 1/2 and 3/2 are the same as in the classical theory of 
Wagner |0J. 

Here is a more detailed account of our calculation of the scaled distribution function. We chose for sampling 13 time 
moments in the interval 120 < t < 2900. The domain statistics is obviously better at earlier times of this interval, 
and it deteriorates at later times, as many domains shrink and disappear. On the other hand, the dynamic exponent 
1/2 shows up, with a good accuracy, only at relatively late times (see Figs. |I] and ^|). Therefore, we had to include 
the relatively late times in our sampling, which led to relatively big error bars in Fig. ^. 

The area fraction e = 0.25, used in our simulations, is moderately large. Therefore, one could expect significant 
deviations of the scaled distribution function, found numerically, from the Wagner distribution |l9| corresponding to 
the same area fraction (that is, having the same second moment). The Wagner distribution has the following form: 

$w(0=Ce ^-^expf (14) 

for £ < v2j and &w(Q = for £ > \/2. The normalization constant 

C = 57 : —- — -T ~ 16.961 , 

Trpe 2 )- 1 +Ei{-2)} 

where Ei{...) is the exponential integral function |32|. 

The two distributions, ^ num {C) ancl ^w{0i are shown in the same Fig. ^. One can see that the difference between 
them is enormous (in order to show the Wagner distribution on the same graph with $„ m , we had to multiply it by 
a factor of 0.5). Therefore, at moderate area fractions, the Wagner's theory is inapplicable. 
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It is instructive to compare the zero moments Mo of the two distributions. The zero moment is the amplitude of 
the scaling law for the number density of domains at large times: n{t) = Mot -1 . We obtained M = 4.72 • 10~ 2 for 
&num ancl Mo ~ 1.43 • 10 _1 for &w- Therefore, for e = 0.25 the Wagner distribution overestimates the number of 
domains at late times by a factor of 3. An additional difference is the pronounced tail in $ num which extends much 
further than the edge of the compact support of the Wagner distribution. Coalescence provides a natural explanation 
to these two facts: coalescence events reduce the total number of domains and produce domains of progressively larger 
size. We shall see in the next Section that an account of coalescence leads to a good quantitative agreement between 
theory and simulations. 



IV. THEORY OF OSTWALD RIPENING WITH COALESCENCE 



In this Section we present a new theory of the globally-conserved (interface-controlled) Ostwald ripening that 
accounts for coalescence. One of our assumptions is that each domain can be represented by an equivalent circular 
domain, or droplet, the area of which is equal to the area of the domain. We shall denote by f(R, t) the distribution 
function of the droplets with respect to their radii. f(R,t) is normalized by the condition J Q f{R 1 t)dR = n(t), 
where n{t) is the number density of the droplets. We start with a brief review of the "classical" theory that neglects 
coalescence and goes back to Wagner [[l9|. Then we derive a kinetic equation that accounts of coalescence. We shall 
focus on the long-time, self-similar asymptotic solutions to that kinetic equation, find the solution by an iteration 
procedure and compare it with the result of the phase-field simulations. 



A. Ostwald ripening without coalescence: a brief review 

At a late stage of coarsening |-ff(t)| *C 1, so there is no nucleation of new domains. Then, neglecting coalescence, 
one can write a simple continuity equation in i?-space for the size distribution function of domains, or droplets: 

d t f + d R (Rf) = Q. (15) 



When criterion ( |5l| ) is satisfied, the dynamics are describable by the area-preserving motion by curvature (|l2j) (where 
we put 5 = 1). This leads immediately to 

R=-^-- --, (16) 
R c (t) R' y ' 

where the time-dependent critical radius R c {t) = \/2/(3 \H (t)\) is determined, at a late stage of coarsening, by the 
conservation of the total area of the minority phase: 



7T / R z fdR = e = const. (17) 



o 



Equations (|15|) -(p^) represent the classical model of interface-controlled Ostwald ripening. This model was formulated 
by Wagner [p_9| by analogy with the theory of Lifshitz and Slyozov (]l8) developed for the locally conserved (diffusion- 
controlled) dynamics. Using Eqs. (|l5|)-(|l7j), one obtains 

R c (t) = J L R !™ = (R(t)) , (18) 

Jo f dR 

where (R(t)) is the time-dependent average radius of the droplets. 

Droplets with R > R c (t) grow at the expense of droplets with R < R c (t) which shrink. The late-time asymptotic 
behavior described by Eqs. (J15|) - (|T7|) is the following [Q. The critical radius grows with time (this corresponds to the 
decrease with time of the effective magnetic field which plays the role of supersaturation) . As a result, a droplet which 
was growing at an early time begins to shrink at a later time. Since all the quantities are position-independent, this 
model represents a mean-field theory. It should be noticed that the mean-field approximation is much more accurate 
for the globally-conserved (interface-controlled) Ostwald ripening than for the locally-conserved (diffusion-controlled) 
Ostwald ripening First, in the globally-conserved case, the "mean field" H(t) is the actual field in the system. This 
is in contrast to the diffusion-controlled Ostwald ripening [|l8| , where a mean field description of the supersaturation 
is an approximation valid only when the typical distances between the droplets are very large compared to the typical 
droplet radius. The second difference concerns the role of correlations. In the locally-conserved case, correlations 
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between droplets result from the Laplacian screening effect, and their relative contribution to the size distribution 
function is of order e 1 / 2 (see, e.g. Ref. The effect of coalescence scales like e (see below) so, at small e, 

correlation effects should be much less significant. By contrast, in the interface-controlled case direct correlations 
between droplets are exponentially small, and significant correlations can be caused only by coalescence events. 
Therefore, in the interface-controlled case, it is legitimate to account for coalescence while neglecting correlations. 

Wagner (lj]] obtained a self-similar solution to Eqs. (|l5|)-(|l7|) (the Wagner distribution) that corresponds to a 
long-time asymptotics of the initial- value problem. The similarity Ansatz is 

W*) = ^*/»(^r)» «■(*) = 

where (3 is a constant number. The scaled distribution ^!f){£) obeys an ordinary differential equation: 



-| + (3 - ij *' (O + (-^ + ± ) * (L> = n . (20! 



The total area conservation ( |l7| ) leads to normalization condition 

POO 

* £, 2 ^p(0d^ = s = const. (21) 
Jo 

Formally solving Eq. (20), one actually obtains a family of solutions parameterized by {3. For a/2 < (3 < 2 a/2/3 



these solutions have compact support: they are positive on an interval < £ < £ m ax(/3), and zero elsewhere. Similar 
solutions in 3D were investigated in Refs. |5|, |2^, p3] |. We call these solutions localized. For < (3 < a/2 the solutions 



of Eq. (20) are extended: they have an infinite tail. These solutions can be written as const ■ x f'o/3(£)i where 



*o/»(0 



(C 2 -2/3^ + 2)2 



exp ( j£ arctan f £ ) . (22) 

Extended solutions fall off like £ -3 as £ — > oo. As a result, the integral in Eq. ( ^l| ) diverges logarithmically, so 
the extended solutions are non-normalizable. Still, as we shall see, they play a crucial role in the theory of Ostwald 
ripening with coalescence. 

Which of the similarity solutions is selected by the dynamics (that is, represents a long-time asymptotics of the 
initial value problem)? It turns out that selection is "weak", that is, determined by the initial conditions. The Wagner 
distribution is selected for (normalizable) extended initial distributions. On the contrary, if the initial distribution 
f(R,t — 0) has compact support, one of the localized distributions is selected. The selection is determined by the 
asymptotics of f(R, t = 0) near the upper edge of its support || |2^, |33|. 

However, this weak selection rule was obtained in the framework of the classical formulation of the problem, Eqs. 
(|H|)-(|T7|). One can expect that strong selection (independent of the initial conditions) can be obtained if one goes 
beyond the classical formulation. Indeed, it was shown in Ref. |3J] (see also [^5[ ) that an account of fluctuations leads 
to strong selection. Fluctuations produce a tail in the time-dependent distribution function and drive the solution 
towards the Wagner distribution. We shall see in the following that an account of coalescence also leads to strong 
selection, even in the absence of fluctuations. 



B. Kinetic equation with coalescence 

We shall now take into account the processes of binary coalescence. Coalescence events occur when two droplets 
contact each other. Within the framework of the GLE, the positions of the droplet centers remain fixed. Therefore, 
for coalescence to happen, at least one of the droplets must be expanding. Consider a droplet of radius R\ < 
R < Ri + ARi. The number density of such droplets is f(R±,t)AR±. Now consider another droplet of radius 
i?2 < R < R2 + Ai?2 in the vicinity of the first droplet. If Ri + R2 > then, during the time interval At, the distance 
between the boundaries of these droplets will decrease by (R± + R,2)At. If the distance r between the centers of the 
droplets obeys the double inequality 



(Ri +R2) <r< (i?i + Ba) + (i?i + R 2 )At , 



(23) 
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(which assumes that the condition Ri + R2 > is fulfilled), then these two droplets will collide during the time 
interval At. Therefore, for the two droplets to collide, the center of the second droplet should be located within a 
circular ring, concentric with the first droplet, with radius R\ + R2 and width (i?i + R 2 )At. The area of this ring is 
equal to 



2tt (R 1 +R 2 )(R 1 +R 2 )At . 
Hence, the average number of such second droplets is equal to 

2M(R 1 ,R 2 )f(R 2 ,t)AR 2 At , 

where 

M(R U R 2 ) =ir(Ri + Ra) {Ri + R2) 0{Ri + R2) 
The total number of the collision events per unit area is equal to 

[2 M{R X , R 2 ) f(R 2 , t)AR 2 At] f(R 1 ,t)AR 1 . 



(24) 



(25) 



Each collision leads to coalescence: disappearance of a droplet of radius R\ and a droplet of radius R 2 , and creation 
of a new droplet. 

Now we make two assumptions that will enable us to construct a closed theory. First, we assume that the area of 
a new droplet, formed by a binary coalescence event, is equal to the sum of the areas of the two merging droplets. 
Second, we assume that new droplet instantaneously becomes circular [ p6[ , so its radius is (R\ + R 2 ) x / 2 . The kinetic 
equation for the size distribution function includes the rates of gains and losses of droplets by coalescence. This leads 
to the following equation: 



dtf + d R (Rf) 



{2M(R u R 2 )x 







8(R - Ri) + 6{R -R 2 )-S[R- J Rf + R 



f(R u t)f(R 2 ,t) alR 1 a!R 2 }, 



(26) 



where S (...) is the Dirac's delta-function and the factor 1/2 is introduced in order to avoid counting each coalescence 
event twice. Performing integration with S(R— Ri) and 5(R—R 2 ) and taking into account the symmetry of M{R\, R 2 ) 
under a transposition of its arguments: 



we obtain: 



M(R u R 2 ) = M(R 2 ,Ri), 



dtf + d R (Rf) 



-2f(R,t) J M(R,R 1 )f(R 1 ,t)dR 1 + 



M(R 1 ,R 2 ) SIR- jRf + R% x 







f(R 1 ,t)f(R 2 ,t)dR l alR 2 . 



(27) 



Integration of the right hand side of Eq. ( j2jj ) over R 2 dR yields zero, so the new kinetic equation preserves the 
conservation law (fT?]) as it should. In addition, the simple relation (|l8|) continues to hold. Integrating the right 
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hand side of Eq. (B7]) over dR, and over RdR : respectively, one can show that the coalescence term reduces the 
number density of the droplets and the total interface length. Moreover, the new equation preserves the dynamic 
scaling. Indeed, if some f(R,t) and R c (t) give a solution to Eqs. (]27|), ( |l7| ) and (Q, then f'(R,t) = rj 3 f(r]R, rft) and 
R' c {t) — r]~ 1 R c (r] 2 t) give another solution to the same equations. This invariance under a stretching transformation 
implies the existence of a self-similar solution that will be considered in the next subsection. 

While deriving Eq. (|27|) , we neglected effects of interactions of three droplets. By this we refer to cases where there 
are three closely lying droplets. In these cases triple coalescence events may occur. In addition, an excluded area 
in the ring (|2^) appears. The effects of this excluded area, and of the triple coalescence events were not taken into 
account in our theory. These effects are expected to be of order e 3 , while the effects of binary coalescence events are 
of order e 2 . Therefore, Eq. (^) is expected to be valid for small area fractions e. We shall see, however, that a very 
good accuracy is obtained even for the moderate value of e = 0.25 used in our simulations, when triple coalescence 
events do occur (see Fig. ||). 

Another limitation of our theory concerns the large- R tail of f(R,t). The tail shape is affected by higher-order 
coalescence events unaccounted for in our theory. This limitation is not very important in practice. The main 
contribution to the critical radius R c comes, for normalizable distributions, by the "body" of the distribution function, 
rather than by the tail. 

We conclude this subsection by a brief discussion of a different type of coalescence: Brownian coalescence. Following 
the pioneering work of Smoluchowski |37| , Binder and Stauffcr ]38| suggested a mean- field scenario of phase separation 
in alloys in which clusters of the minority phase are regarded as Brownian particles: they perform random walk in 
space. When two clusters collide, they merge into a larger single cluster. The corresponding kinetic equation includes 
an integral term whose general structure resembles that of the integral term in Eq. (^7|), but with a different kernel 
M(Ri, R2). If the cluster diffusivity is a power-law function of the cluster size, one arrives at a self-similar solution 
for the size distribution function of the droplets. An important further development was the work of Siggia J3S| ] 
who considered hydrodynamic interactions between randomly moving and coalescing droplets in phase separating 
binary fluids. Following the work of Siggia, the Brownian coalescence in binary fluids has been extensively studied 
theoretically and experimentally. Among important issues here is a crossover from Ostwald ripening (the Lifshitz- 
Slyozov- Wagner mechanism) to Brownian coalescence jj], |39[ , plethora of hydrodynamic interactions in the process 
of coalescence |}9| ^l], [42^1 , scaling violations jl3) etc. In parallel, Brownian coalescence has been investigated in the 
context of coarsening of clusters of atoms or vacancies diffusing on surfaces, following particle deposition [Q. It is 
clear that Brownian coalescence is different in its nature from the coalescence process considered in this work. In 
contrast to Brownian coalescence, droplets in our system do not move: they coalesce only because they grow. 

C. Self-similar solution with coalescence 
Equations (^7|) and ( |l7| ) admit the same similarity Ansatz as Eqs. ( |l5| ) and (p7|): 

/(*•*) = ^*(S?) ( 28 ) 

and 

R c (t) = p- 1 1 1 / 2 , (29) 

where /3 is again an unknown yet constant number. The scaled distribution function <&(£) obeys the following integro- 
differential equation: 



-I + "-£H {)+ (-H 



-2*(0 J w(t,ti)*(£i)d£i+ 







^i>6)«U-VS + 6 d6 (30) 
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subject to normalization condition 



In Eq. (p0|) we denoted 



/>OC 

Jo 



(31) 



2/3 



1 1 



a 6/ ' 

where 9 (...) is the theta-function. It is convenient to rewrite Eq. ( |30| ) in a symbolic form: 

£/3<I> = Afp[$\ , 

where 



■| + /3-^(0+(-| + ^]«f.i.. 



(32) 



(33) 



(34) 



and 



[$] (o = -2$(o y tflfofosfcodei- 









^(a,e 2 )5^-^i 2 +e 2 2 ) 



(35) 



One important property of Eq. Q3C|) can be noticed immediately: the coalescence term vanishes identically at 
< £ < 1/(2/3). As a result, the scaled distribution function at < £ < 1/(2(3) should coincide (up to a £- 
independent multiplier) with one of the solutions of the classical Wagner's problem. A simple argument shows that 
parameter 0. parameterizing this solution, should be less than \/2. Indeed, inverting the linear operator Cp, we 
rewrite Eq. ( p3| ) as an integral (rather than integro-differential) equation: 



*(0 = *p(0 



(36) 



where functions &p(£) were introduced in subsection A and C\ is a constant. Unless < y/2, the integral over d£' 
diverges. Therefore, *£p(£) should be one of the extended solutions ^0/3 gi yen by Eq. (f22|). In addition, since the 
second term in the square brackets of Eq. ( |36| ) would lead to divergence of the integral appearing in the normalization 
condition ([jl|), we must choose C\ — 0. Hence, Eq. ( ^6| ) reads: 



*(0 = <M0 



(37) 



Integral equation (|3^ ) and normalization condition ( pi] ) make a complete set. For a given e, the scaled distribution 
function $ = &p(£,) and parameter — 0(e) are uniquely determined. Therefore, an account of coalescence does 
provide strong selection to the problem of Ostwald ripening. 
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D. Solving Eqs. (37) and (31) 



Our procedure for solving Eqs. ( |37| ) and (31) employs the one-to-one correspondence between e and /3. Therefore, 
one can fix (3 and solve Eq. ([37]) by iterations for <f> = $ l a(^). The normalization condition d3l]) is not used at this 
stage. After a sufficiently accurate estimate for $ = <&[}(£) is obtained, one employs Eq. (RW to calculate e that 
corresponds to this f3. Repeating this procedure for different [3, one obtains the family of solutions ^(C) and the 
dependence e = e(/3). Inverting this dependence, one arrives at (3 = /3(e) and finds the correspondence between e and 
scaled distributions <&/?(£) = $p( E )(0- 

Now we introduce an iteration scheme that implements this idea. The scheme exploits the fact that, for the exact 
solution of Eq. @, 



*(0=X*o/j(0 at 0<^< 



1 

2/3' 



(38) 



where \ — x(P) is ^-independent and unknown in advance. Let us introduce an auxiliary unknown function 0(£) 
X^^iO !§■ By definition, 



HO = *o/3(0 at < £ < 



In terms of the new function <j>, Eq. (p7h becomes 



2/3 



(39) 



, v |-/3 + i)^(a 



The iteration scheme for Eq. (fftT) is the following: 



(40) 



(|-/3 + ^)* 0/3 (e)' 



(41) 



for A; = 0, 1, 2, . . . . We use the arbitrariness of \k and demand that, at each iteration, Eq. (B9|) is satisfied: 



&+i(0 = *o/j(0 at 0<^<^ 



(42) 



This implies 



Xfe 



1/(2/3) 



(|-/3 + i) * 0/3 (0 



(43) 



for k = 0, 1, 2, . . . . Equation ( ff3| ) and ( fll| ) define the iteration process explicitly. If the sequence <pk ,k = 0,1,2,... 
converges to a finite limit 0(£), then the sequence Xk , k — 0, 1, 2, . . . converges to a finite (positive) number Xi and we 
can find $(£) = X<?K£)- What is left is to use Eq. ( |3l| ) and compute the corresponding e. If the convergence of the 
iteration scheme is fast enough, then 



®k(0=XkM0 



(44) 



and 



£fc(/3) = 7r Xfe / eMOdH, fc = l,2, 



(45) 



give a good approximation to the solution already after a small number of iterations. 

We implemented this iteration procedure numerically. As it is clear from Eqs. (41 ), (E3) and (|45|), the numerics 
involve only calculations of (double and triple) definite integrals. We started with the trial function 0o(O = ^0/3 (0- 
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The advantage of this trial function is that it already satisfies Eq. (|39|). We performed a detailed investigation of 
the convergence of this scheme for different values of (3 up to iteration k = 8. The results of this investigation (and 
proofs of convergence of the integrals at infinity) will be presented in a separate publication. The main result is that, 
with this choice of the trial function, the convergence is very fast in the body of the scaled distribution function. For 
example, in the case of e = 0.25 the first iteration already gives an accurate result. Convergence in the tail is a more 
subtle issue that will be addressed separately. 

The /3(e)-dependence, found by this procedure is shown in Fig. [j]. One can see that, as e — > 0, [3 — > \[2 from below. 
Overall, this behavior is expected. Coalescence effects become small at small area fractions and, as e — > 0, the scaled 
distribution function should approach the Wagner distribution. It is surprising, however, that for quite small e (for 
example, 0.01), the value of (3 is still significantly different from V2- 

Figure H shows the scaled distribution function for (3 = 0.93 (which corresponds to e ~ 0.25). The solid line shows 
the result of the first iteration < &i(f)- This result agrees well with the phase field simulations, even at the moderately 
large value of e — 0.25. This is a strong evidence in favor of the major role of coalescence in the interface-controlled 
Ostwald ripening. The dashed line is the trial function < E'o(^)- One can see that this trial function (with (3 — 0.93) 
already gives a fairly accurate estimate to the solution in the body of the scaled distribution. Therefore, the non- 
normalizable extended distributions ( |22| ) do play a special role in the theor y. In the tail region the solution falls off 
more rapidly than \&o^t so there is no problem with normalization condition (pit). The zero moment Mq corresponding 
to the once-iterated numerical solution (the solid line in Fig. [|) is equal to 5.08 • 10~ 2 which is in good agreement 
with the value 4.72 • 10~ 2 obtained in the simulations. 

An independent estimate of (3 at e = 0.25 is provided by the dynamics of the effective "magnetic field" H(t). The 
area- preserving dynamics, Eq. (|l2|), imply that |£f(f)| — |\/2k/3|. Assuming that all droplets have circular shape 
(an assumption already used in our theory), we have |k| = l/(R(t)). Employing the similarity solution, we arrive at 
|k| = 13/t 1 / 2 . Therefore, l/\H(t)\ ~ ct 1 / 2 , where c = 3/(V2f3). For the corrected power-law fit shown in Fig. || c = 2.3. 
This corresponds to (3 ~ 0.92 which is remarkably close to the value of 0.93 obtained by the iteration procedure. 

V. CONCLUSIONS 

We investigated globally conserved phase ordering dynamics in systems with short range correlations. The numerical 
simulations were done with a 2D phase-field model: a Ginzburg-Landau equation with a global conservation law. The 
sharp-interface limit of this equation: the area-preserving motion by curvature was introduced, and a criterion for its 
validity formulated. Assuming dynamic scaling within the model of area-preserving motion by curvature, we obtained 
the 1/2 dynamic exponent for the characteristic coarsening length scale (for critical and off-critical quench), and 
for the effective "magnetic field" (for off-critical quench). Our numerical simulations for critical quench and for an 
off-critical quench with area fraction e — 0.25 confirm these scaling laws. The results for critical quench coincide with 
earlier results, obtained by particle simulations. The scaled form of the equal-time pair-correlation function is found 
to weakly depend on the area fraction, similar to the locally-conserved systems. 

Recently, dynamic scale invariance and "normal" scaling have been reported in the same globally conserved coarsen- 
ing system, but for long-range (power-law) correlations in the initial conditions 0. Our present results, combined with 
those of Ref. [Q, indicate "normal" scaling properties of globally conserved systems for any generic initial conditions. 

The main focus of this work was on the dynamics of the size distribution function of the minority phase domains 
for the off-critical quench. Our phase field simulations show that this distribution function exhibits, at large times, 
a self-similar form. The scaled distribution function, however, is dramatically different from the well-known Wagner 
distribution, despite the fact that correlations between domains are negligible in globally-conserved systems. We 
attributed this strong deviation to coalescence that provides, at small e, a leading correction to the Wagner's theory. 
We developed a new theory of Ostwald ripening that takes into account binary coalescence events. The theory 
possesses dynamic scale invariance and yields a non-linear integro-differential equation for the scaled distribution 
function. For a given area fraction, the problem has a unique solution. Therefore, coalescence renders a strong 
selection rule to the problem of Ostwald ripening: the scaled distribution is selected by the area fraction of the 
minority domains. The scaled distribution function predicted by the new theory is in good agreement with the scaled 
distribution obtained from the phase field simulations for a moderate area fraction of 0.25. In addition, the new 
theory accurately predicts the amplitude of the late-time power law of the effective magnetic field Hit). Deviations 
from the classical Wagner's theory remain significant even for very low area fractions. Therefore, coalescence plays a 
major role in the interface-controlled Ostwald ripening. 

We hope that the results of this work will stimulate further experiments on globally-conserved interface-controlled 
systems. 
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Appendix. Criteria for the area-preserving dynamics 



Let us find the general conditions under which the area-preserving dynamics [Eq. (|l2|)1 represent an accurate 
approximation to the more general sharp- interface theory [Eqs. @ and (R)]. We first notice that in the general 
sharp-interface theory it is the quantity A(t) — L 2 H(t)/4, rather than A{t), that is conserved ||. Hence, a necessary 
condition for the validity of the area-preserving dynamics is simply 



H(t) < 



L 2 



(46) 



Additional criteria are found in the following manner. We notice that Eq. (Q) includes the same combination 
Sk+(Sii) 1/2 gH(t) as the one that appears on the right hand side of Eq. (||). Therefore a constancy or, more precisely, 
slow variation of A(t) implies a slow variation of H(t). Correspondingly, the term on the left hand side of Eq. (J7]) 
should be small in this case compared to each of the two terms on the right hand side. We can exploit this fact and 
formally solve Eq. (Q) perturbatively: H(t) — H^(t) + h(t) where H (0 \t) = -g^ 1 (S/ n) 1 ^ 2 k, is the leading term, 
and h(t) is the subleading term. Substituting Hit) into Eq. (|7|) and keeping terms up to order h(t) we obtain: 



h(t) = 



L 2 



4g(6n) 1/2 A(t) 



£T(°)(t) = - 



VM(t) 



k(s, t) . 



(47) 



The perturbation expansion is valid as long as 



\h(t)\ < |J? (0) (*)I and 



\h(t)\ 



(6fj,y/ 2 A(t) 



\h(t)\ 
L 2 



Using the expressions for H^{t) and h{t) we see that the dynamics described by 
area-preserving if the following two inequalities hold: 



and (0) are (approximately) 



n(s,t) 



(<$Ai) 1/2 A(t) 



< 



and 



KM)I 



L 2 



(48) 



d_ 

dt 



k{s,t) 



(Sfj,) 1 / 2 A(t) 



k(s,t) 



(49) 



Now we can see why the area-preserving dynamics do not apply to a single circular droplet. In this case the zero 
order of the perturbation theory would give k(s, t) = 1/R = const, where R is the droplet radius. Therefore, criterion 
(||) is violated. 

Therefore, neither of the two criteria 



For a "donut" one has n(s,t) = 0. Therefore, neither of the two criteria ( [48| ) and (|49j) is obeyed, and the dynamics 
are not area-preserving. If the system consists of Ni(t) domains of the majority phase and N2 (t) domains of the 
minority phase, the Gauss-Bonet theorem |4q] yields 



k(s, t) 



MjViCt)- N 2 (t)) 

A(f) 



(50) 



where N\(t) does not include the large "sea" of the majority phase. We shall employ this relation in the following. 

Let us check criteria (filf), ( [48| ) and ( |49| ) for the "standard problem" of phase-ordering, when the initial conditions 
describe a disordered state with short-ranged correlations. We assume that the system exhibits, at late times, dynamic 
scale invariance. In other words, the characteristic domain size l(t) is the only relevant length scale. For an off-critical 
quench, the phase-ordering morphology is that of Ostwald ripening: competition of droplets of the minority phase, see 
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below. This implies the following scaling relations: n(s,t) ~ l/l(t), k(s,t) ~ 8/l 3 (t), and A(t ) ~ (£)/(£) ~ eL 2 /lit). 
In addition, if(t) ~ A/Z(t) and A = eL 2 . Using these relations in any of the inequalities j4q), fl48| ) or (49), we arrive 
at 



A 



< e. 



(51) 



That is, the more general sharp- interface theory is reducible to the area-preserving motion by curvature as long as 
the ratio between the interface width and the coarsening length scale is much smaller than the area fraction of the 
minority phase. If the coarsening system remains a two-phase system, and if the system size is big enough, this 
condition holds, at late times, for any nonzero area fraction 47 . Notice that, at e <C 1 it takes more time for the 
system to reach the area-preserving regime. 



A. J. Bray, Adv. Phys. 43, 357 (1994). 

L. Schimansky-Geier, C. Ziilicke and E. Scholl, Z. Phys. B-Cond. Mat. 84, 433 (1991). 
J. Rubinstein and P. Sternberg, IMA J. Appl. Math. 48, 249 (1992). 

C. Sire and S.N. Majumdar, Phys. Rev. Lett. 74, 4321 (1995); Phys. Rev. E 52, 244 (1995). 

B. Meerson and P.V. Sasorov, Phys. Rev. E 53, 3491 (1996). 
A.D. Rutenberg, Phys. Rev. E 54, 972 (1996). 

A. Peleg, M. Conti, and B. Meerson, Phys. Rev. E 64, 036127 (2001). 

P. Wynblatt and N.A. Gjostein, in: Progress in Solid State Chemistry, edited by J.O. McCaldin and G. Somojrai, (Perga- 
mon, Oxford, 1975) Vol. 9, p. 21. 

M. Zinke-Allmang, L.C. Feldman, and M.H. Grabow, Surface Sci. Reports 16, 277 (1992). 
Z. Toroczkai and E. Williams, Phys. Today 52 (12), 24 (1999). 

I. S. Aranson, D. Blair D, V.A. Kalatsky, G.W. Crabtree, W.K. Kwok, V.M. Vinokur, and U. Welp, Phys. Rev. Lett. 84, 
3306 (2000). 

I. Aranson, B. Meerson, P.V. Sasorov and V.M. Vinokur, e-print cond-mat/0107443 

A. J. Bray, Phys. Rev. Lett. 66, 2048 (1991), and references therein. 

J.F. Annett and J.R. Banavar, Phys. Rev. Lett. 68, 2941 (1992); L.L. Moseley, P.W. Gibbs, and N. Jan, J. Stat. Phys. 
67, 813 ( 1992). 

B. P Lee and A.D. Rutenberg, Phys. Rev. Lett. 79, 4842 (1997). 
J.A. Glazier, M.P Anderson, and G.S. Grest, Philos. Mag. B 62, 615 (1990). 
G.S. Grest, M.P. Anderson, and D.J. Srolovitz, Phys. Rev. B 38, 4752 (1988). 

I.M. Lifshitz and V.V. Slyozov, Zh. Eksp. Teor. Fiz. 35, 479 (1958) [Sov. Phys. JETP 8, 331 (1959)]; J. Phys. Chem. 
Solids, 19, 35 (1961). 

C. Wagner, Z. Elektrochem. 65, 581 (1961). 

I. Aranson, B. Meerson, and P.V. Sasorov, Phys. Rev. E 52, 948 (1995). 
M. Marder, Phys. Rev. A 36, 858 (1987). 

A.S. Mikhailov, Foundations of Synergetics I. Distributed Active Systems (Springer- Verlag, Berlin, 1993). 

A. Peleg, B. Meerson, A. Vilenkin, and M. Conti, Phys. Rev. E 63, 066101 (2001). 
M. Gage, Contemp. Math. 51, 51 (1986). 
S.M. Allen and J.W. Cahn, Acta Metall. 27, 1085 (1979). 

B. Meerson and I. Mitkov, Phys. Rev. E 54, 4644 (1996). 
W. Ostwald, Z. Phys. Chem. 34, 495 (1900). 
T.M. Rogers and R.C. Desai, Phys. Rev. B 39, 11956 (1989). 

P. Fratzl, J.L. Lebowitz, O. Penrose, and J. Amar, Phys. Rev. B 44, 4794 (1991). 
K. Humayun and A.J. Bray, J. Phys. A 24, 1915 (1991); Phys. Rev. B 46 10594 (1992). 
K.R. Elder and R.C. Desai, Phys. Rev. B 40, 243 (1989). 

S. Wolfram, The Mathematica Book, 3rd ed. (Cambridge Univ. Press, 1996), p. 1081. 
B. Giron, B. Meerson, and P.V. Sasorov, Phys. Rev. E 58, 4213 (1998). 
B. Meerson, Phys. Rev. E 60, 3072 (1999). 
I. Rubinstein and B. Zaltzman, Phys. Rev. E 61, 709 (2000). 

In fact, new droplets just formed by coalescence do not have a circular form (see Fig. ^). However, for small e, deviations 
from the circular form are significant for big droplets, that is in the tail of the size distribution function. In the "body" of 
the size distribution function (which gives the dominant contribution to the critical radius R c ), the droplets are already 
close to circular. 
[37] M. von Smoluchowski, Z. Phys. Chem. 92, 129 (1917). 
[38] K. Binder and D. Stauffer, Phys. Rev. Lett. 33, 1006 (1974). 
[39] E.D. Siggia, Phys. Rev. A 20, 595 (1979). 
[40] C.K. Haas, J.M. Torkelson, Phys. Rev. E 55, 3191 (1997). 



15 



[41] H. Tanaka, J. Chem. Phys. 105, 10099 (1996). 

[42] V.S. Nikolaev, D. Beysens and P. Guenoun, Phys. Rev. Lett. 76, 3144 (1996). 
[43] A.J. Wagner and J.M. Yeomans, Phys. Rev. Lett. 80, 1429 (1998). 

[44] C.R. Stoldt, C.J. Jenks, P.A. Thiel, A.M. Cadilhe and J.W. Evans, J. Chem. Phys. Ill, 5157 (1999); D. Kandel, Phys. 

Rev. Lett. 79, 4238 (1997). 
[45] One can check that a straightforward iteration scheme for would diverge. 
[46] E. Kreyszig, Differential Geometry (University of Toronto Press, Toronto, 1959). 

[47] A criterion similar to Eq. ( ^l| ) also appears when the initial condition represents a single long stripe [^3). In that case l(t) 
is the stripe length. 



16 




FIG. 2: Snapshots of globally conserved coarsening for off-critical quench with area fraction e = 0.25. The upper row corresponds 
to t = 5.2 (left) and 32.3 (right), the lower row to t = 204.8 (left) and 1305.5 (right). 
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x=r/l(t) 



FIG. 3: Scaled correlation function obtained by numerical simulations with the GLE subject to a global conservation law 
for times t > 15. The solid line is C(x) for critical quench, and the dashed line is C(x) for off-critical quench with area 
fraction e = 0.25. The dotted (dashed-dotted) line represents C(x) obtained in particle simulations of globally conserved || 
(non-conserved J3FJ) dynamics for critical quench. 
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FIG. 4: Characteristic coarsening length I vs. time for the critical quench (circles), and off-critical quench with area fraction 
e = 0.25 (squares). The two solid lines are corrected power-law fits l{t) = lo + bt a with a — 0.50, lo = 0.5, and b = 1.2 for the 
critical quench, and a = 0.51, lo = 1.3, and b = 0.9 for the off-critical quench. The dotted line shows a pure t 1 ^ 2 power law. 
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FIG. 5: l/\H(t)\ vs. time for an off-critical quench with area fraction e = 0.25 (circles). The solid line is a corrected power law 
fit l/\H(t)\ =a + ct a with a = 7.4, c = 2.3, and a = 0.51. The dotted line is 2.3 1 1/2 power law. 
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FIG. 6: Scaled distribution function of domain sizes <3>(£), where £ = R/t 1 ^ 2 . The diamonds represent $ num (^): the scaled 
distribution obtained in the simulations with the GLE for times t > 120. The error bars show (twice) the variance of the scaled 
distribution functions for 13 time moments in the interval 120 < t < 2900. The dotted line is the Wagner distribution &w(0> 
Eq. ([lij), for the same area fraction e = 0.25. In order to show it on the same graph with <& num we had to multiply it by 0.5. 
The dashed and solid lines show the distributions $o(f) and respectively, predicted by the theory of Ostwald ripening 



with coalescence, 
[see Sec. 0, Eq. (I 



IV 



resented in Sec 

for fi — 0.93 which corresponds to e 



These distributions represent the zero and first iterations of the iteration procedure 
0.25. 
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FIG. 7: Parameter j3 versus e (solid lines) as predicted by the theory of Ostwald ripening with coalescence. This dependence 
was obtained by a single iteration applied to Eqs. ( |37| ) and (^l|). The dashed line shows the limiting value j3 = %/2 expected at 
e — > 0. The solid circle is the point found in our phase-field simulations. Figure (a) shows £ in a logarithmic scale, figure (b) 
in a linear scale. 



